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Abstract 

On the basis of previous work by Fadin, Lipatov, and collaborators, and of our 
group, we extract the "irreducible" part of the next-to-leading (NL) BFKL kernel, 
we compute its (IR finite) eigenvalue function, and we discuss its implications for 
small-x structure functions. We find consistent running coupling effects and sizeable 
NL corrections to the Pomeron intercept and to the gluon anomalous dimension. 
The qualitative effect of such corrections is to smooth out the small-x rise of struc- 
ture functions at low values of Q^. A more quantitative analysis will be possible 
after the extraction of some additional, energy-scale dependent contributions to the 
kernel, which are not treated here. 
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The small-x rise of structure functions at HERA |l|] has stimulated an impressive 
theoretical effort [2-15] in order to understand the high-energy behaviour in QCD. 

In a series of papers [2-3,5-9,11], Fadin, Lipatov and collaborators have investigated 
the high-energy cluster expansion (Fig. 1) of the parton-parton cross-section, with the 
purpose of generalizing the leading-logs BFKL equation[0] to the next-to-leading (NL) 
order. Similar results for the (heavy) qq production cross section]^] and eigenvalues |]^, [T2[], 
and for the squared gluon emission amplitude [|10| have been produced by other authors. 



The outcome of such an analysis is the calculation of some " irreducible" vertices which, 
defined by gluon Regge-pole factorization [Q] , have the role of incorporating the low energy 
features of the QCD scattering amplitudes. They are, more precisely, the two-loop gluon 
trajectory renormalization [§, the one-loop reggeon-reggeon-gluon (RRG) vertex [5-7] 



and the RRQQ § and RRGGR K^ clusters at tree-level (Fig. 1 (d-f)). 



Such vertices, which involve parton counting in the final state, suffer from mass sin- 
gularities and need be combined in a sum, with parton number n = 0, 1,2 , in order to 
define the IR finite, irreducible part of the NL kernel. This cancellation of singularities 
was shown in Ref. [ISj . 



Recently, the A^j-dependent (or qq) part of the kernel was extracted by the authors 
|1T2| , |l3,its eigenvalue was computed, and the ensuing Pomeron shift and anomalous di- 
mensions were evaluated. Both shift and resummation effects turn out to be suppressed, 
in this case, by a nonplanar colour factor, which needs not be present in the gluonic case. 

Interest in extracting the full gluonic contribution is thus substantial, but is not as 
simple as in the qq case, however. In fact, in order to define the irreducible vertices in the 
cluster expansion, we need to subtract the leading kernel iteration with logs accuracy. 
In addition, we need to define an off-shell scale for the energy, starting from a partonic 
cross-section which is not an IR safe hard process and thus suffers from Coulomb-like and 
possibly coUinear singularities due to initial partons. 

In performing this procedure, we shall distinguish two kinds of contributions to the 
kernel: (i) the properly irreducible ones, depicted in Fig. 1, which come from the RR, 
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RRG and RRGG clusters to be defined below, and (ii) some additional contributions, 
that we call energy-scale dependent, which are remainders of the leading term, with its 
energy scale, after the required factorization of Coulomb and coUinear singularities. The 
last step involves the choice of a factorization scheme of vertices and kernel which should 
allow the use of the latter in hard processes. 

The purpose of the present note is to perform in part this program, by combining 
the "irreducible" terms mentioned above and by discussing their eigenvalues and related 
features. This will allow us to understand the running coupling effects and the main 
consequences for high energies and anomalous dimension behaviour. On the other hand, 
full quantitative results can be obtained only after the extraction of the energy scale 
dependent terms, which is deferred to a subsequent analysis. 

Let us start by defining the irreducible terms more precisely. We work in cij-space 
(Mellin transform in the energy variable s) and transverse momentum space with respect 
to the incoming partons' axis. The leading kernel in D = 4 + 2e dimensions has the 
well-known form 

oj ujL [1 — e) [ki — k2j^ uj 

(1) 

where we have adopted the notation of Fig. 1 (i. e., k's (q's) for the exchanged (emitted) 
momenta), and we have introduced the one- loop gluon trajectory p, 



and the notation 



a, = ^, „,(/<^) = ^^r(l-.). (3) 



We also understand that transverse integrations carry the measure (i[k] = ci^^^+^^k/ 



TT 



The kernel (1) is the prototype of the e-dependent kernels to be written out. It is 
finite for e and q = ki — k2 7^ 0, but may still have singular eigenvalues because 
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of the q = singularity. The virtual term - in this case the one-loop gluon trajectory - 
regularizes the singularity by providing a subtraction which yields a finite eigenvalue. 

In fact by applying the kernel (4) to the test function (kl)"'^^, with < 7 < 1, and by 
using the integral 



d[k2 



r(l -e)q2 Vkf 



2\ 7-1 



kfy i r(i + 5)r(7 + e)r(i-7-e) 

e r(l-e)r(7 + 2e)r(l-7) 



e 



expiexoil)) + :,e\^P'{l - 7) - 3^'(7)) + 0{e') 



it is easy to combine real and virtual terms to obtain the characteristic function 



Xo(7) = 2^(l)-^(7)-^(l-7), 



(4) 



(5) 



where iplj) is the logarithmic derivative of the F-function. 

The NL kernel contains virtual and real emission terms also. The virtual term is the 
two-loop gluon trajectory 

2cu(2)(k2^ 



K 



(NL) 



-TT 



UJ 



(6) 



which, by collecting the gluonic contributions only, is given by 



4e2 



-2 / 3 



TC 



6 



67\ /202 IItt^ ,^ . 



1^ 



(7) 



^/iV V 6 V 6 18/ V 27 

The subtraction of the leading term (~ 2{uj^^^ log s)^) is in this case unambiguous, because 
the only available scale is k^. 

The real emission terms involve several scales, and thus some prescription is required 
for the subtraction of the leading term. We start by writing the one-gluon emission 
amplitude with the one-loop corrections of Ref. [5-7] as follows 

MiI) = Mi°) fl + /3jHk?)+/3«(k^) + 

S2 



-^(^)(k?) f log ^ + log + -^(^)(k^) flog 4 + log \ + Mi^) 



(8) 



where si, S2 are subenergy variables of the emitted gluon of polarization e and colour 
a, (3^^\(^B^) one-loop corrections p, ^ to the A (B) vertex, M^^^ is the leading 

amplitude, and M*-^-* is defined to be the NL irreducible one. They are given by 

s 



k?k: 



PA • q 



pb • q 



(9) 



and by 

Re MiP 

k? 



T 



{'^9,e ■ J) 



C(e) TTcosTre / q' 



sm Tie 



.2\ £ 



11 



k2 
^2 



+ 



3 (k? - ki)2 V 6 k? 



Ilk? + k2 



k^kg 



k^ 3 (k^ k2)^ 



loi 



k? 

ki. 



+ 



k^k^ 



3(k? 



ki) 



+:^(i-(2q^-k? 



11 + q- 



k? + k2 



2q2-k?-k^ 



(k? - k 



2N2 



k2 



where J'^ is the gluon-emission current associated with high energy scattering 
Jf* = — ^ is the soft insertion current, occurring in M^^^ only. 



(10) 
while 



PB-q 



Notice that we have used in Eq. (10) the small-q^ behaviour of Eq. (21) of Ref. 
together with the fixed q^ form of Eq. (86) of Ref. [Q. Furthermore, we have incorporated 
the u^^^ logk^ terms in the definition of the leading part in Eq. (8), so as to subtract them 
out in a scale invariant form. 

The one-gluon contribution to the NL kernel comes from the interference of M^^^ in 
Eq. (10) with M^^^ in Eq. (9). Using the polarization sums 



-J' 



4kfk| _ 4ki ■ k2 



q^ q^ 

and performing, for simplicity, an azimuthal average of the polynomial part in q^, we 
obtain 



K 



(NL) 
l9 



uj 4 



C{e) {q^/fi'^y IT cos ire 
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r(i-e) 
1 



q- 



£r(i-£)q2 kf-k 



sm ne 

.2 



2 u2l0gT-2 

2 



1 



k^ 

2 

2 / J 



(12) 



Finally, the two-gluon emission cluster has been recently computed in Ref [|TT|, where 
the authors suggest subtracting the leading term 

d[qi] f dx 



(13) 



q^(q — qi)2y x{l — x) 

with a scale invariant rapidity phase space J x{i^x) ~ 21og(l/5). By using this prescrip- 
tion, we obtain, from Eq. (20) of Ref. the expression 



K. 



(NL) 
29 



u 4 



2\e 



[ r(i-e)q2 

-ifeoZz(kl,k2) +if(ki,k2 
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¥ ^ lis 



TT 



/202 



7C(3) e 



+ 



(14) 



where we have introduced the "collinear" kernel 



Hcoii(ki,k2) = — 



K2 K]^ 



X 



1 



log tan — + Jm Li2 I i 
v/kfk! V k^ 



|ki| 



|ki| 



vr 

314^ 



-2 ^ 

in which an azimuthal average has been performed, and the dilogarithmic one 



^(ki,k2) + 



vr 



2q ■ (ki + ks) 



3k2 



+Li2 1 
»i 



k? 



q2(ki + k2)2 
Lt2 i 1 



^2 

log ^ lo, 



k?k 



21^2 



k2 



Li 9 



2 (k? + ki)2' 
k? 



k2 , 

*^2/ J 



+ 



+ 2 

+(ki < — ^ 



dt 



tk,^2 



L^o (ki - t«^2 

Ml- 



q^ 



q2(ki 



k,)2 



t)k 



X 



kl{l-t) + qH 



(15) 



(16) 



Let us now investigate the physical features emerging from the irreducible NL kernel, 
as defined by the sum of Eqs. (6), (12) and (14). In order to find its eigenvalue, we 
shall proceed as for the leading term by applying the kernel to test functions of the form 
(kf)^"^, with < 7 < 1. Due to its explicit renormahzation scale dependence we expect 
the outcome to contain factors of log(k^//x^). 

We combine first the e-dependent singular terms in Eqs. (12) and (14) and we notice 

that the most singular ones (~ 1 / s'^) cancel out directlyj^ leaving the total real-emission 

^This justifies performing the expansion up to relative order e^, instead of e^. 
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singular part 



11 



ri-e) 



11 

~6 



18 ^ "6 



202 
"27 



7C(3) e'+ 



+ 



3 eq^T{l-e)_ 

This kernel has a finite e ^ limit at fixed q^, but its eigenvalues are still singular. 



(17) 



We then compute the eigenvalues of the kernel (17) by using Eq. (4) and the additional 
integral 



ri-e) 



2\ 7-1 



klY' 1 r^(i + £)r(7 + £)r(i-7-2£) 



£2 r2(i -£)r(7 + 3e)r(i -7) 



kfV' 1 



exp (2exo(7)) + tt^' (4V^'(1 - 7) - 8^'(7)) + 0{e^) 



and, by combining them with the virtual term (6) we obtain the finite result 



u 4 



xo(7) ( -y f + Y - y) - y (^o(7) + x;(7)) + 6C(3) 



(19) 



Note now that the coefficient of the log /i term is precisely the leading kernel eigenvalue 
with a beta-function coefficient. Therefore, it can be interpreted as a running coupling 
factor, much as for the qq contribution. We can thus express the total (L+NL) kernel in 
the form 



1 - 6a,(/i2) log 4 ) iro(ki, ka) + a,(/i')iri(ki, ks 
/i 



(iro(ki,k2) + a.Ki(ki,k2)) 



(20) 



which defines the NL scale-invariant kernel K^. 



Factorizing the running coupling at the scale is an asymmetrical procedure, but 
is convenient for the discussion of the non-scale-invariant BFKL equation [|13|. Using a 
different scale (e. g., as(k?,)) implies changing K\ so as to leave the total NL kernel 
invariant. 

Finally, a straightforward calculation allows the computation of the characteristic 
function of the remaining finite part of the kernel, except for iJ, whose eigenvalue is 
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estimated semi-analytically to be (7) — Yln=i^n[{l + n) ^ + {l — '-f + n) ^],withai = .72, 
a2 = .28, 03 = .16. We thus obtain the gluonic part of the Ki eigenvalue in the form 



y (xo(7) + Xo(7)) + (y - y ) ^0(7)+ 



+ 6C(3) + 



TT 



37(1 - 7) 



TT 



cos7r7 



sin 7r7 / 3(1 — 27) 



11 + 



7(1-7) 



;i + 27)(3-27) 



(21) 



to be compared with the iV/-dependent part obtained previously |T2 



1 5 1 

o W(7) + Xo(7)) - 0X0 - ( - 



7r 



cos 



7r7 1 + 17(1-7) 



3"" A/'2Vsin7r7y 1 - 27 (1 + 27)(3 - 27) 

(22) 



Correspondingly, the (azimuthal averaged) k-space kernel can be rewritten in a more 
compact form: 



(9) 
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- H,,u{k,, ka) + ^(ki, ks) + 67rC(3)5(2)(q) 



R 



(23) 



where we have used the notation of Eqs. (15) and (16), and we have introduced the 
regularized distributions in 2-dimensional transverse space 

/(ki,q)|^ = /(ki,q)e(q2-A2)-52(q) / f{k„q)d\ (24) 

JA2 



The parallel expression for the gg-part, obtained in ref [0, has the form[| 



(</) _ Nfas 



6tt 



5\ 1 



- ^i/afc(ki, ks) 



where the abelian contribution is defined by 



Hab{ki,k2) = — 



(25) 



k2 

+ (22-?i-?| 



k?, -i|k2 



log —2 tan 



k2 



|ki 



+ Im Li2 ( i 



|k2| 
|ki| 



(26) 



^Eq. (3.17) of Ref. contains some misprints, in particular a C/ factor instead of Nf. 
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We notice that in both cases the natural scale of as appears to be q^, rather than , and 
that the "coUinear" and "abelian" terms have same sign and similar structure, but that 
the latter is suppressed by the l/N^ colour factor. 

Several comments on the results in Eqs. (20)-(26) are in order. Consider first using 
them to describe the gluon density in a two-scale hard process. According to our general 



arguments |jT3|, |Tj] the kernel with running coupling is consistent with the renormalization 



group for kf > k2 > A"'. 

More precisely the Mellin transform Gi^{Q,Qq) of the BFKL Green's function is rep- 
resented by 



Gu>{Q,Qq) = (exp 

7+V-X (7+) V 

in the anomalous dimension regime 



^ -f+{as{t'))dt'^K{u,to) (27) 



t = log ^ > to, bujt>x[^], (2^ 



A2 ^' ^ V 2 

where 7+ ~ 7gg + §j7gg is the larger eigenvalue of the anomalous dimension matrix, 
defined at NL level by 

l = ^(Xo(7+) + «.Xi(7+)). (29) 
(jj 

From the definition (29) both perturbative and resummed expressions of the NL 
anomalous dimension follow from the formula 

Xi (7l (ff)) 

where 72,(0:5/0;) = as/^ + 0{as/uj)'^ is the well-known leading gluon anomalous di- 
mension. 

Therefore, the low order expansion of 7+ 



7+ = - + aJAi + A2— + A3 — +... (31) 

UJ \ OJ \UJ / J 

implies the small-7 behaviour of Xi 

Xi(7)^^ + — + ^3 + 0(7), (32) 
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which can be checked on Eqs. (21) and (22) to be consistent with the known [|l^ expres- 
sions in the DIS scheme 

llN^as Njas 1 . _ Nja, fh 13 



127r 67r N^' 67r V3 6A/^?y 

For the gluonic part, Hcou is responsible for the colhnear behaviour, because H vanishes 
for either k2 = or = 0. 

Furthermore, the gluonic eigenvalue provides, through Eq. (30), important resumma- 
tion effects which are driven by the negative value A^^^ = — llAi'c/127r at the 7 = 1 double 
pole in Eq. (21). This term causes a rapid increase of ^— 7^Xi^''(7) j relative to its value 
at 7 = (Fig. 2(a)), which reaches a factor of 2.5 at 7 = 1/2. Thus, unlike the qq part, 
the gluonic part is expected to be relevant for scaling violations at HERA. 

The second important point concerns the high-energy behaviour expected for the gluon 
density. We pointed out in Ref. Jl^ that in the running case, two kinds of critical 
cu-values occur. One is the singularity of the anomalous dimension expansion occurring 
in Eq. (29) close to 7 = 1/2, for which we get an ag-dependent value 



= ctsXo (1 - aas) . (34) 

The other is the true Pomeron, the w-singularity dominating the high-energy behaviour 
beyond the anomalous dimension regime. While the former admits the rough estimate 
(34), the latter turns out to be dependent on the behaviour of close to k2 = A^, and 
thus cannot be really predicted. 

If we take the formula (34) as a qualitative estimate, we realize that the NL gluon con- 
tributions in Eq. (22) yield a rather large negative shift, namely a ~ 3.4 with our present 
knowledge (Fig. 2(b)). This would mean that the "Pomeron" intercept is substantially 
decreased, of the order of ^ .2 for = .15. 



Let us emphasize that this indication cannot be taken yet as a quantitative estimate, 
because of the scale-dependent contributions to the kernel that we have neglected and 
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because various cross-checks of the whole approach are still needed. It means, however, 
that the NL corrections go in the direction of bridging the gap with soft physics, by 
smoothing out the small-x rise at low values of Q^. 

If the above magnitude of NL corrections is confirmed, it raises the problem of the 
slow convergence of resummed perturbation theory at small-x. Fortunately some classes 
of corrections can be roughly understood at all orders, because they correspond to physical 
phenomena we already know about. 

One class of corrections is due to the collinear behaviour of large-x contributions, 
which in the cluster expansion approach occur in higher order clusters, and give rise to 
multiple poles of xil) at 7 = and 7 = 1. Resumming these poles is mandatory ||T^ to 
understand the lower eigenvalue of the anomalous dimension, and in general the behaviour 
of xin) close to 7 = 0. To see the point, we have plotted in Fig. 4 the function 

2^(1) - ^(7 - Aias) - ^(1 - 7 - Aias) + a. \xi{l) - M f-^l ) (35) 

\ \sm 7r7 / / 

which coincides up to NL level with xq + Q^sXi) but differs at higher orders, by a resum- 
mation of the collinear behaviour. 

It is apparent that the effect of resummation is to displace the 7 = singularity to 
7 = A\as < 0, as expected, and to reduce the Pomeron shift, by about 15% for as = .15. 



Another class of corrections was noticed long ago by one of us ||T^ to be due to co- 
herence effects in the soft gluon emission region. The ensuing structure function equation 
with angular ordering, further investigated by Catani, Fiorani and Marchesini [jl8| (the 



CCFM equation), is the basis for the treatment of such effects to all orders. 

From the results which are already available [|l^, it appears that coherence plays a 



role starting from the constant A3 in Eq. (32), and thus from the three loop level. It 
constitutes, therefore, a rather delicate check of the whole approach, because the scale- 
dependent terms, neglected here, are also expected to contribute. 

On the whole, we think that a full understanding of the next-to-leading kernel will 
put several phenomenological issues on quantitative grounds and will help to bridge the 
gap with large-x properties, low physics, and diffractive phenomena. 
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Figure captions: 



Figure 1: (a-c) Leading and (d-f) next-to-leading contributions to the high energy 
cluster expansion. The loop number of radiative corrections is specified, and clusters with 
external particles are omitted. Intermediate quarks and gluons are understood in (f). 

Figure 2: Plot of (a) the gluonic contribution to 7^X1(7) (b) the cts-dependent 
" Pomeron" . 

Figure 3: Plot of the resummcd characteristic function in Eq. (35), with the sym- 
metrical choice Xiil)^ fo^^ which Q;s(k>) is factorized. The corresponding Pomeron shift 
should be decreased by Aa;p = — lasCtsTr^ to compare with the result in Eq. (34). 
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